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1 Introduction 
1.1 Motivation 

This paper is motivated by the problem of source localization using a large sensor network. In this context, the 
observation is a complex valued M-variate time series (yn)nez represents the number of sensors of the array) 
given by 

K 

Yn = L ^k.n^^djc) + V„ = A(0)s„ + V„, 
fc=l 

where 
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• The K < M scalar (in general complex valued) time series ((Sjt „) „^z for k- l,...,K are non observable, and 
represent the signals transmitted by K transmitters. The vector s„ is given by - . . . , SK,n) ^ ■ 

• For each k, 6]c is a scalar real parameter characterizing the direction of arrival of transmitter k. 9 ^ a(0) 
is a known C*^-valued function depending on the sensor network geometry, and matrix A{0) is defined as 
A(0) = (a(0i),...,a(0if)). 

• (Vn)nEZ finally represents an additive complex Gaussian noise, i.e. v„ - (fi, yM,?;)^ where the M time se- 
ries {{V]c^n}n€z) ic=i m mutually independent identically distributed (i.i.d.) sequences such that Re(ffc,„) 
and lm{vk,n) are independent real Gaussian random variables with zero mean and variance a^/Z. 

The classical source localization problem consists in estimating vector - {6\,...,9k)^ from N samples collected 
in the M x N matrix Y^f - (yi,...,yiv). This problem was extensively studied in the past (see e.g. |21| and the 
references therein). The so-called subspace estimator ofO- {9\,...,9k)^ is based on the observation that if ma- 
trices A(0) and Siv - (si, . . have both full rank K, then the angles {9k)k=\,...,K are solutions^ of the equation 
a(0)*njva(0) - 0, where 11^ represents the orthogonal projection matrix on the kernel of matrix A(0)SArS^A(0)*. 
The existing subspace methods consist in estimating for each 6 the quadratic form rj^id) - a(0) * njva(0) of Hn by 
a certain term f]Ni9), and then to estimate the K angles as the argument of the K most significant local minima of 
function 9 — ► f\Ni9). This approach has been extensively developed when N — ► +oo and M fixed. In this context, 
riN{9) can be estimated consistently for each by Tyjv (0) = a(0) * njva(0) with n jv the orthogonal projection matrix 
on the eigenspace associated to the M - K smallest eigenvalues of the empirical covariance matrix -^YivY^. It 
clearly holds that supg^^.^j^i \fiNid) -ilNi6)\ converges forwards almost surely, and this allows to prove that the 
corresponding estimators (0fc)fc=i,...,js: of the direction of arrivals are consistent. 

If however M and N are of the same order of magnitude, a quite common situation if the number of sensors M 
is large, then the above estimators show poor performances because IIjv is no longer an accurate estimator of IIjv. 
In order to study this context, Mestre & Lagunas [18] were the first to propose consistent estimators of 7]^ (0) when 
M,N —- +00 in such a way that cjv = — ► c, with c > 0. In Mestre & Lagunas [18], it is assumed that the source 
signals {Sk,n)k=i,...,K are mutually independent complex Gaussian i.i.d. time series with unit variance elements. 
Under this assumption, y^, can be written as 

where Ry - A(0)A(0) * +o^\m represents the covariance matrix of the time series (y„) „ez and x„ is a complex stan- 
dard Gaussian vector. Matrix 11 ]v coincides with the orthogonal projection matrix over the eigenspace of Ry associ- 
ated to the eigenvalue o^, and Mestre & Lagunas addressed the problem of estimating consistently any quadratic 
form of n^v from the empirical covariance matrix -i-YjvY^ where Yjv - Ry'^X^v and Xjv - (xi, . . . ,xjv). Mestre & 
Lagunas 1 1 8 1 used properties (see Silverstein & Choi 1 20 1 , Bai & Silverstein |T1 f2l ) of the empirical covariance ma- 
trix, and were able to exhibit a M x M matrix tlnd.N such that a^ftj-j-rf jv^at - a^IIjvajv for each deterministic 
bounded sequence of vectors (ajv) when M,N ^ +oo in such a way that cjv = ^ c, with c > 0. In some sense, 
matrix n,-,,^ TV can be viewed as a consistent estimate of IIjv but in a weak sense because in general, it does not hold 
that ||n,-,rf TV ~ IIjv II ^ 0, where we have denoted by ||.|| the usual spectral norm. Mestre 8d Lagunas concluded that 
for each 9, a(0)*n,-,rf Tva(0) is a consistent estimate of rjNid)- However, the consistency of the angular estimates 
was not established. Note that these results do not require any hypothesis on K which may scale with N or not. 

In Vallet et al. L23J, a more general case was considered where the time series (sj; „)fc=i ^ are deterministic 

signals for which the spectral norm of matrix -^A(0)Sjv is bounded w.r.t. the dimensions M,N,K. This time, 
random matrix Yjv is non zero mean, and corresponds to the so-called "Information plus Noise model" investi- 
gated in various works of Girko (see e.g. 1 1 1 1 ) and Dozier & Silverstein ( 1 9 1 , 1 8 1) . Using new results on the almost 
sure localization of the eigenvalues of the empirical covariance matrix -^YjvY^, |23| generalized the estimator of 
Mestre & Lagunas 1 18 1, and derived a "weakly consistent estimator" 11^ of IIjv, i.e. a^nwajv - a^IIjvajv for 
each deterministic bounded sequences of vectors (ajv). Therefore, it holds that for each 9, fjjv(0) - a(0)*njva(0) is 
a consistent estimate of rjN {,9) if a(0) is uniformly bounded in N. 

The goal of the present paper is to pursue the work |23|, and to establish that the angle estimates defined as 
the K most significant local minima of function 9 a(0) * flNaid) are consistent. As it will be shown below, the 

^ The K angles are the unique solutions under certain assumptions on function 6 — » a(0) 
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consistency of the angle estimates is based on the property 

sup \f]N{9)-VNm\^0 (1) 

d£[-n,n] 

almost surely, that we shall refer to as the uniform consistency of the estimate fjNiO) of rjNiO). 

This paper is organized as follows. In Section|2] we provide some background material on the asymptotic eigen- 
value distribution of the large information plus noise model, on the almost sure localization of the eigenvalues of 
the empirical covariance matrix, and on the consistent estimator of a*^IlNaN proposed in (23). In Section|3l we 
prove the property of uniform consistency of estimator 77 jv (0) (see (1)) when function a(0) is defined by 

a(0).^fl,e'0,...,e'f^-i'«]^ (2) 

(1) of course holds for more general functions, but we believe that considering the typical example defined by (2) 
is informative enough. The proof of (1) heavily relies on results concerning the probability that the eigenvalues 

of ^'^^ escape from the intervals in which they are located almost surely for N large enough. These results are 
believed to be of independent interest. Finally, we establish in Section|4]the consistency of the K most significant 
local minima of function Q fiN{6) by following the approach in fTSl. 



1 .2 General notations and useful results 

We now introduce various notations and results used throughout the paper. 

• If £ c IR, Int(ij) and dE represent the interior and the boundary of E respectively. 

• If z £ C, the complex conjugate of z is denoted z or z*. For a complex matrix A, we denote its transpose by 

and its Hermitian adjoint by A* . 

• We denote by °°(IR,[R) (respectively S^'^(IR,[R)) the set of all smooth real-valued functions (resp. compactly 
supported smooth real values functions). 

• The quantity C will represent a generic positive constant whose main feature is to be deterministic and 
independent of M and N. The value of C may change from one line to another. 

• Similarly, Pi and P2 will denote generic polynomials, independent of M and N, with positive coefficients. 
The polynomials may change from one line to another. 

• Complex Gaussian distribution: A complex valued random variable Z - X + iY follows the distribution 
'^jV [a + ip, fj^) if X and Y are independent real Gaussian random variables jY{a, ^ j and ^ j re- 
spectively. The variance of Z, denoted as Var(Z) is defined as Var(Z) -E\Z -E[Z]\^ - . 

• Poincare inequality (see Chen \ G\y. let Zi - X\ + iY\, Zp - Xp + iYp be p iid '^JY{<d,a^) random variables 
and consider a function 7 defined on IR^^ continuously differentiable with polynomially bounded partial 
derivatives. Then, if X = {Xi,... ,Xp)^ and Y = (Yi, . . . , Yp)^, the random variable 7(X,Y) can be written as 
7(X,Y) = 7(Z,Z) and 

_ P 

Var [7(X,Y)] = Var [f{Z,Z)] < (J^ ^ 

where we define as usual the differential operators ^ = | [ ^ ~ ' ^ j and ^ = |[^ + '^]-IfTis real- valued, 

it is clear that '^^i-'^ coincides with the complex conjugate of ^^^"^ . In this case, the Poincare inequality 
reduces to 

Var(7(X,Y)) <2cr^^E 



dfiZ.Z) 



dzi 



df{Z,Z) 



dzi 



df{Z,Z) 
dZi 
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Stieltjes transform: Let be a positive finite measure on K. Its Stieltjes transform m is the function defined 
by 

m{z]- \ , Vz £ C\supp(fi), 

Jk A - z 

where supp(/i) represents the support of measure /x. Function m is holomorphic on C\supp(/i) and satisfies 
> for z £ C\IR and m{iy) when y ^ +00. Moreover, supp(/z) c 1R+ if and only if ^'^[^^f^ > for 
z £ CW. The mass of the measure /i can be evaluated through the formula 

ii(IR) = lim -iym{iy). 

y— *+oo 

We also notice that if m(z) is the Stieltjes transform of positive measure /i, then it holds that 

\m{z]\ < ^ r < 



and that m'(z) - satisfies 



dist(z,supp(/x)) |Im(z)r 

. catiofioc 

a-z)2 

\miz)\< ^< 



dist(z,supp(fx)) llm(2)l 
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on C\supp(/i). We finally recall the following version of the inverse Stieltjes transform formula: For each 
function y/ e ^~(IR,IR)), we have 

f i//(A)d^(A)= -limlmff y/{A)m{A+ iy)dA]. (3) 
JR yio \JR I 

2 Background on the Information plus Noise model and on the estimator of 

f23] 

All along this paper, we consider integers M,N,KeN* such that 1 < < M, = K[N) and M - M{N) are functions 
of N with cn - ^ ^ c> as N ^ 00. We assume that 

Assumption A- 1 : < cat < 1 andO < c < 1. 

In this section, Z]v represents the complex valued M x N random matrix given by 

2jv=-;= = Bjv+Wiv, 

where Bjv = ^^^=^ and W]v = Matrices Bjv and Wjv are assumed to satisfy the following assumptions 
Assumption A-2: Matrix Bn is deterministic and satisfies sup^^ || Bjv || < +00 
Assumption A-3: Rank(B jvB^) - K < M where K may scale with N or not 

2 

Assumption A-4: The entries ofmatrixWff are i.i.d and follow the complex normal distribution '^jY{Q, 

We assume moreover that the non zero eigenvalues of B]vB^ have multiplicities 1 in order to simplify the 
notations. In the following, we denote by = Ai,jv - Am-k.n < Am-jst+i.at < ■■■ < Ajvf.jv and (Ujt,jv)fc=i,...,M the 
ordered eigenvalues and associated eigenvectors of BjvB^. The eigenvalues and the eigenvectors of matrix Zjv2)^ 
are denoted (A^ ^v)fc=i,...,M and (Ufc,Af)fc=i,...,M. and fiN represents the empirical eigenvalue distribution of Zjv^^^ 
defined by 

Y M 
P-N- — y ^} 

As we assume cn < I, the joint probability distribution of (Ajt jv)*:=i,...,m is absolutely continuous (see e.g. James 
[16)) and it holds that the (Aj; ^v) fc=i,...,M have multiplicity 1 almost surely. We finally denote by Qjv(z) the resolvent 
ofmatrixZjvX^, i.e. Qjv(z) = (ZjvZ^ - zIm)"^ 
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2. 1 The asymptotic eigenvalue distribution /i at of jUjv 

It is well-known ([TT] Th. 7.4], (5] Th. 1.1]) that it exists a sequence of deterministic probability measures (/ijv) such 
that fiN - pjv — weakly almost surely. Measure hn is characterized by its Stieltjes transform rriNiz) which is 
known to satisfy the equation 



1 

mjv(z) = — Tr 



-z(l + a'^Ct^mj^[z])lM + 0-2(1 - Civ)lM + 



1 + a^CNniNiz) 



(4) 



for each z e C\[R^ . In the following, we denote by S^n the support of /ijv. As fiN - fiN —'nO weakly almost surely, it 
holds that 

mpjiz) - mpjiz) (5) 
almost surely for each z £ C \ IR^ . The following result wiU be of help. 

Lemma 1 ([HI, (4]). Lety/ e Sf^(lR, K) and (r^f) a sequence of holomorphic functions on C\U such that 

|rjv(z)|<Pi(|z|)P2[--4^1, 
lllm(z)| j 

with Pi and Pz two polynomials with positive coefficients, independent ofN. Then, 



limsup / i//(x)r]v(x+ iy)dx <C<oo, 

ylO MR 

with C a constant independent ofN. 

Taking into account the previous result, it is shown in |23 1 that 

E[m]v(z)] = mjv(z) + -^^^, (6) 

with as in Lemma [1] Using the inverse Stieltjes transform formula @. we obtain that for each function i// e 
'^?°{IR,IR), it holds that 



1 1 *^ » r II 

— E[Tri/.(2iv2:;v)] = ]^ZE[vUfc,Jv)]= n>{X)pN{dX) + 0\j^ 



(7) 



If we denote by T]v(z) the matrix- valued function defined by 
Tiv(2) = 



z(l -KT^C]vmjv(z))lM + cr-'d - Cjv)Im + ■ 



1 -I- a^CNniNiz) 



then Tat coincides with the Stieltjes transform of a positive matrix valued measure /i^ with support 5^jv such that 
fi^iS^N) = Im (see Hachem et al [HI Th. 2.4 & Prop. 2.2]), i.e. 



T.(z).[ 



As mjv(z) verifies the equation (4), it is clear that j^Jr^^ - Pn- 

In the remainder of the paper, we will make use of the following result proved in (23] if Wjv is complex Gaussian 
and in Hachem et al 1 14 1 in the non Gaussian case. 

Theorem 1. Consider two sequences of deterministic vectors (b]v), (djv) such thatsupp^ llbjvll < +oo andsupj^ Wd^W < 
+00. Then, it holds that 

b;vQjv(z)djv-b;rTiv(z)djv — 0, (8) 

N 

almost surely for each z e 



5 



2.2 The characterization of the support of /iat 

The support S^jq of was first studied in Dozier & Silverstein |8| and a more convenient cliaracterization was 
presented in [23]. We first recall (see |8|) that if z £ converges forwards xeU, then, m]v(z) converges torwards 
a finite limit stiU denoted mjv(x). Function x rriNix) is continuous on U, continuously differentiable on U\dS^N, 
and verifies Eq. (4) on R\dS^M. Moreover, /Xjv is absolutely continuous and its density coincides with function 
ilm(mjv(x)). 

In order to present the chacterization of S^n, we first introduce the following notations. We denote by /iv,0jv 
and the functions defined by 

/Af(w)= ■^Tr(BjvB^- wIm)~^, 

(pNiw) = w(l-0-^Cjv/jv((A'))^ + (T^(l-Cjv)(l-0-^Cjv/jv(w)), 

wjvU) - z[\ + a^Cf]mN{z)f- - a'^[\ - Cf])[l + a^Cf]m]^{z)). (9) 
We are now in position to characterize 5V. 

Theorem 2. The function cj)N admits 2Q non-negative local extrema counting multiplicities ( with 1 < Q < K + I) 
whose preimages are denoted W^^ < < w^^^ < w^j^j-- - £ i^qm < n- Define - (pNiw~j^) and jc^^ - 
</>jv(if^^) forq=l...Q. Then, 

■"■l.N ^ -"-l.N - -"-Z.N < ^Q,JV < ^Q,N' 

and the support S^N of p^ is given by 

Moreover, for q - 1, . . . , Q, each interval 1 jv' '^(/ iv ^ contains at least an element of the set {0, Am-k+\,n, ■ ■ ■ > ^m.jvI 
and each eigenvalue of B^B*^ belongs to one of these intervals. 

The second statement of the theorem shows that each eigenvalue of Ba/B^ corresponds to a certain interval of 
S^N- More precisely, an eigenvalue of BjvB^ will be said to be associated to cluster [x~ j^,x'^ ^] if it belongs to the 
interval {w~ j^, w'^ j^)- We note that the eigenvalue is necessarily associated to the first cluster [x^ j^,x^ j^]. 

We finally recall the useful properties of function w^ defined by (9) (see L23J). We still denote by w^ix) the 
limit of w^iz) when zeC^ converges torwards xeU. 

Proposition 1. Function wn : C ^ C satisfies the following properties: 

• Function x W]\iix) is continuous on U. and continuously differentiable on U\d5^N, 

• Im(ifAr(z)) > (flm(z) > 0, 

• W]\i is real and strictly increasing on IR\^]v, 

• WNix-j^) = w~j^ and wn{x'^j^) = w+ ^ for each \<q<Q, 

• Im(ifAr(x)) > if and only ifx e Int(5^]v). 

2.3 Some useful evaluations 

In this paragraph, we gather some useful bounds related to certain Stieltjes transforms. We first recall that the 
inequality 

|1 + a^CNmN(z)\ > Red + a'^ c^mNiz)) > 1/2 (10) 

holds for z £ C (see Loubaton & VaUet [13). We now consider function z — — Proposition 2.2 in 1 13| 

implies that it coincides with the Stieltjes transform of a probability measure carried by K^. Moreover, fTO) shows 
that the support of this measure is included in 5^jv U {0}. Therefore, we obtain that 

1 |z| 

5 < (11) 

|l + CT'^C]vmjv(z))| |Im(z)| 
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for each z e C\R as well as 

1 \z\ 



|l + o-2cjvmAf(z))| dist(z,5^Af) 

1 that r 

We now claim that the inequality 



for each z £ C*\5V- We also recall that matrix TAf(z) satisfies TAr(z)Tjv(z)* < -^-^ for z e C+ (see 1 13, Prop. 5.1]). 



Tiv(z)TAf(z)*< (12) 
dist(z,.5>)'^ 

also holds on C\5^jv- In order to establish {12), we follow the proof of Proposition 5.1 in J3). We first remark that 
fimction miv(z) defined by 

1 - Cjv 

mjv(z) = CNiriNiz) 

z 

is the Stieltjes transform of probability measure jl^ - cni^n + (1 - cjv)(5o. The support of /ijv thus coincides with 
5^N U {0}, and is included in IR+ . Therefore, it holds that '""i^u)'^" > if z £ C\B«. We remark that 

T]v(z)-T]v(z)* , , , f dfij^W 
Im(z) / 



2i 'J^^,|A-z|2' 
By using the identity, Tjv(z) -Tiv(z)* = Tjv(z) (Tjv(z)"* -Tiv(z)"i)TAf(z)*, we get after some algebra 



Im(z) 



^„ |A-zH 

|1 + (t2c^^jy(z)|2 



Im(z)Tiv(z)Tiv(z)* +(j2lm(zmA,(z))TA,(z)TA,(z)* + ^^—^ ^Im(miv(z))Tjv(z)BjvB;vTjv(z)*, 



for each z £ C\1R, or equivalently 

f diij^m ^ Im(zmAf(z)) Q-^Civ Im(mAf(z)) 

/ =Tjv(z)TAf(z) +0- — — Tjv(z)Tjv(z) +- ^ — ^ — — — — Tjv(z)BjvB^Tjv(z) . 

Jsf^ |A-zr Im(z) 1 + (t2c^^^(z)|^ \m[z] 



\X-z\ 

Consequently, we obtain that 

rfMjv(A) 



Tjv(z)Ta,(z) 



Js 



.y^ |A-z|2 



for z £ C\IR, but also for z £ C\5% because both members of above inequality are continuous on C\5V- This 
immedialely leads to (12). This inequality also implies that for each z £ C\^jv, 

min lAfcjv- Wiv(z)| > -dist(z,,S^iv)- (13) 
k=\,...,M ' ' 2 

Indeed, Tjv(z) can be written as T]v(z) = (1 + ct2cjy?7i^(z)) (BjvB^- wjv(z)Im)~^- Therefore, ||Tjv(z)|| is equal to 

NT r.m II + g^Civmjv(z) I 

l|Tiv(z)|| = ^ 1- — T 

mmfc=i,...,M|Afc,iv- wjv(z)| 

so that (13) follows from (12) and (TO) . 

Since mAr(z) is the Stieltjes transform of the distribution ■j^Z^j'^^^^. it holds that 

1 

|mjv(z)| : 



dist(z, {Ai,iv,...,AM,Af}) 
as well as 

1 

Wj^{Z)\ < ; J. 

dist(z,{Ai,jv,---,AM,jv}) 
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We now consider the rational function z ^ - — — — - which will play an important role in the following. Its poles 

are solutions of the equation 1 + a^CNinMiz) - 0, and satisfy some useful properties. From now on, we denote by 
Ajv the diagonal matrix Ajv - Diag (Ai,jv, . . . , ^m,n] and hy SIn the matrix 

njv = Ajv+— -^11^, (14) 
M 

where 1 denotes vector 1 - (1, 1, . . . , 1)^. We denote ft)i,Ar < . . . < cjm.n its eigenvalues. Then we have the following 
straighforward properties. 

• The zeros of z 1 + o^CNniN^z) are included in the set {wi,jv.- ■ ■<'j^m,n}- 

• If the eigenvalues Ai,Af,..., Ajvf,w of 'Ln'L*j^ have multiplicity one, the equation 1 + a^CNtfiMiz) = has M 
multiplicity one solutions which coincide with the {(i)k,N)k=i,...,M- Moreover, Ai w < wi,Ar < ... < Xm,n < 

&M,N- 

• If the eigenvalue Ajt jv has multiplicity p> I, i.e. Afc_i,Af < Aj;^ - Aj;+p_i^ < ^k+p,Ny then, 

^k-l,N < ^k,N - (^k,N - Xk+p-2,N - (>>k+p-2,N - ''^k+p-l.N < ^k+p-\,N < ^k+p,N, 

and the Wj-.iv that do not coincide with some eigenvalues of Zjv2^ are zeros of 1 + a^CNrhNiz). 

Remark 1. Since cjv < 1, we recall that the eigenvalues {^k,N)k=i,...,M have multiplicity 1 almost surely. However, in 
subsection \3.2\ it will be necessary to define properly the solutions of I + a^CNinNiz) - everywhere. This explains 
why the case where some of the iAk,N) k=i,...,M are multiple has to be considered. 

Function z ^ zii+ah-NifiNiz)) Stieltjes transform of a probability measure whose support coincides with the 

set of all roots of the equation z(l + cr^cjvmjv(z)) - 0, which is included into the set {O.wi^jv, . . ., wjvf.iv}- Therefore, it 
holds that 

1 , 1^1 



|1 + a^CNmN{z)\ dist(z, {0,&)i,jv, ■ ■ ■,(^m,n}) 

for zeC\{0,cbi N,---,ii)M n) and 

1 Izl 



ll + CT-^Civmjvtz)! |Im(z)| 

for z e C\IR. We eventually notice that 

1 



(15) 



IIQjv(z)|i 



dist(z,{Ai ,iv,...,AM,Jv}) 



2.4 Almost sure localization of the eigenvalues (Afc^v) k=i,...,M 
We recall the two following useful results of 1 23 1 and 1 1 7 1 . 

Theorems (|23|). Assume assumptions A- 1 toA-4 hold. Let a, beU,e>Q and No e N such that 

]a-e,b + e[nS^N - 0, 

for each N > Nq. Then, with probability one, no eigenvalue of'Lf^l.*^^ belongs to [a, b] forN large enough. 

Theorem 4 (| 17|). Assume assumptions A- 1 toA-4 hold. Let a, b e U, e > 0, Nq e !M such that]a-e,b + e[nS^N - 
for N > Nq. Then, with probability 1, 

cardjfc : Xk,N < a] - cardjfc : Ak,N < w^ia)] (16) 
card{fc : Xk,N > b] - cardjfc : Xk,N > w^ib)] (17) 

forN large enough. 

It is useful to mention that sup^Xg^ < +oo and that these two theorems are still valid ifb- +oo (see 1 171). 
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2.5 The consistent estimate of quadratic forms of Iljv 

Let IIjv be the orthogonal projection matrix on the kernel of BjvB^ and let [an) jven be a sequence of deterministic 
M-dimensional vectors such that sup^ l|aiv|| < oo. Then, |23 1 proposed a consistent estimate of 77 jv defined by 

T]N = a^njvaAT. 

The approach of |23 1 is valid under the following assumptions. 

Assumption A-5: For N large enough, none of the strictly positive eigenvalues o/B jvB^ is associated to the first 
cluster [x^ j^,x^ j^], i.e. Xm-k+i,n > t^Nix^pj] forN largeenough. 



Assumption A-6: It holds that 



< liminf Xj < limsup xj^jy < liminf jv- 

JV — hoo 



iV— '+00 ' jY— »+oo ' A^— *+oo 

Using theorems|3]and|4l we deduce that if t^,t^,t^, t^ are real numbers independent of N satisfying 

< rj~ < liminf xj'jY < limsup xj^jy < < < liminf Xj ^ < limsup xi jv ^ ^2^ ^^^^ 

iV— +00 ' JV— +00 ' JV— +00 ' AJ— +00 

then, almost surely, for N large enough, it holds that 

< fj < Ai_Ar < ... < Ajvf-jc.jv <ti <t2 < Xm-k+\,n < ■■■< ^m,n < ^2- (19) 

Assumptions|5]and[6]thus imply that, almost surely, the smallest M- K eigenvalues of Z jvZ^ are separated from the 
K greatest ones for N large enough in the sense that the 2 sets of eigenvalues are included into 2 disjoint intervals 
that do not depend on N. It is interesting to remark that Assumptions [5] and [6] are "deterministic conditions" 
depending only on a^,CN - ^.and on the eigenvalues of BjvB^. If K remains fixed, recent results of Benaych- 
Rao |3| (see also |17|) imply that Assumptions A-5 and A-6 hold if and only if liminfjv— +00 /ljvf-*:+i,JV > f^^/c. If 
however K scales with N, the derivation of more explicit conditions equivalent to Assumptions [5]and[6]is still an 
open problem. 

We are now in position to present the consistent estimator of tjn proposed in |23] . It is based on the observation 
that 

njv = [ (BjvB;^ - ^Im) " ^ dA, 

where represents a contour enclosing and not the strictly positive eigenvalues of BjvB^, and the symbol 
means that the contour is oriented clockwise. The estimator of |23] is based on the observation that under As- 
sumptions |5] and [6] function ifjv(z) provides such a contour for N large enough. In the following, for y > and 
e > 0, £• < I small enough, we consider the rectangle defined by 

S^y-{z-x+iv,0<t^-3e<x<t^+3e<t2-3e,-y<v<y} (20) 

and its boundary dS^y. Then, the properties of function if jv(z) (see Proposition[T) imply that for N large enough, 
the set w^idSSy) is a contour enclosing the origin, but not the other eigenvalues of BjvB^. Therefore, IIjv can also 
be written as 

njv=— / (BjvB^- wjv(z)Ijw)~^ w^(z)dz 
2i7r Ja%" 

or equivalendy 

1 r u>'{z) 

njv = — / Tjv(z) 2^ — dz (21) 



because (BjvB^ - Wf^{z]\M] ^ - i+^CNmN(z) ' Using l|5) and (8) as well as the following lemma 

Lemma 2. Almost surely, for N large enough, the M solutions (Wfc.jv) k=i,...,M of the equation 1 + a^CNmN(z) - 
satisfy 

h < Ai,jv < (Wl.JV < ... < Ajvf_jf,jv < Ci)M-K,N <t\ <h < ^M-K+1,N < ^M-K+l,N < ■■■< Ajvf,jv < WjVf.JV < ^2 
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It is showed in \23i that matrix YIn defined by 



Un^ — [ 

2l7T Jd 



w'{z) 

QNiz)- 2^^— — -dz (22) 

where wjvC-z) - z{l + a^CNrfiMiz))^ - a^{l - cjv)(l + a^ct^rht^iz)), satisfies a^njvajv - a^IIjvajv almost surely. 
We note that the poles of the integrand of the righthandside (r.h.s) of (22) coincide with the set {Ajt jv.'&fc.iv '■ k - 
l,...,M], which by (19) and Lemma|2l verifies 

dist[dSSy,{Uk,N,<^k,N)k=i,...,M^) > 3c (23) 

almost surely for N large enough. In pratice, the above estimator is quite easy to implement because, as the local- 
ization of the poles of the integrand in (22) w.r.t. the contour d^y is known (see lemma|2), the contour integral in 
(22) can be solved, and expressed in closed form in terms of the (u^^jv, ^k.Ny^^k.N) k=\,...,M- 



3 Statement and proof of the uniform consistency of estimate f) at (0) 

From now on, we assume that vector a(0) is given by (2) and that assumptions[5]and[6]hold. We consider fj", fj^ , 
and satisfying fT8) as well a rectangle S&y defined by (20) . We prove here the following result. 

Theorems. Assume assumptions A- 1 toA-6hold. Then, we have 

sup |f?jv(0)-7?jv(0)| '0. 

ee[-7t,7rl 

with probability one. 

In order to prove theorem [S] we show that it is sufficient to establish that for each a > and for each 6 e 
[-71, n], lP(|a(0)* (fiiv - niv)a(0)| > a) decreases fast enough forwards 0. For this, a tempting choice is to use the 
Markov inequality, and to establish that the moments of a(0) * (fljv - IIjv} a(0) decrease fast enough. However, the 
observation that (23) holds for N greater than a random integer does not necessarily imply the existence of the 
moments of a{9)*YlNaid). In order to solve this technical problem, we establish that the probability that at least 
one element of {Afc jv.'i'fc.iv '■ k- l,...,M] escapes from [t^ - 2e, t^ + 2e] u - 2c, t2 + 2c] decreases at rate for 
any I eN, and prove that the moments of a convenient regularized version of a(0)* (fijv - nAr)a(0) converge fast 
enough forwards 0. 

In the following, we denote by 3~c the set 

- [ti - e, t^ + £]IJ [t2 - £, t2 + e], 

We first establish in Sections |3 . 1 I and l3^ that the events <Si^n and iS2,n defined by 

Si,N - {at least one of the (Aj; ^v)fc=i,...,M escapes from 3^], (24) 
Sz.N - {at least one of the {cbk,N)k=i,...,M escapes from 31]. (25) 

verify P(<S'i,Af) - ^(777) for each / e Using this result, we introduce in Section [3?3l the regularization term, de- 
noted xn, defined as follows. We consider a function <p e 'W^{R, M.^) satisfying 

f 1 for A £ 

(pW^i , ^ ^ . (26) 

[0 forA£R\([r]"-2e,r+-H2e]u[f2 -2e,r+-H2c]) 

and (p{X) e (0, 1) elsewhere, and define the random variable 

XAT = det (I){1.n^*n] det (n jv) , (27) 

which verifies l^c^ < xn where Sn - (oi.jv u S2,n- We will prove that, considered as a function of the real and imagi- 
nary part of the entries of Wjv, Xn is a '^^ function, and using Poincare inequality, we will establish that 

1 



E|a(0)*(niv-njv)a(0)xjv|^'=( 



for each integer /. The above mentioned properties eventually allow to prove the uniform consistency of estimator 

flN{9). 
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3.1 Evaluation of the escape probability of (Afc,iv) k=i,...,N 

The purpose of this section is to prove the following technical result. 
Proposition 2. Under assumptions A- 1 -A-6, for each I e N, it holds that 

P(^i,iv) = ^f-^ 



To prove this result, we consider a function y/o e IR^) such that 

f 1 for A £ 

VoW^-i (28) 
[0 for Ae [fj",r+]U[f2.^2']- 

and y/o{A] e (0, 1) elsewhere. From this definition, we clearly have 

PtA.jv) < P(Tn//o(2jvZ^) s 1) < E[(Tn//o(Zjv2:^))^'] 

for Z E N. In order to establish Proposition[2] it is therefore sufficient to prove that E |(Tri/'o(2iv2;^))^'j - [j^^ 
for each integer I which is the object of the next lemma. 

Lemma 3. Assume assumptions A- 1 toA-6hold. Then, for all function y/ e {U,U) constant over the complemen- 
tary of a compact interval and which vanishes on the support of for all N large enough, it holds that 

E[(Tri//(Z^,Z^)f (29) 

for each / e M. 

Proof: We prove Lemma|3]by induction on /. We first consider the case / = 1, and consider a function iff as above, 
and denote by C the constant value taken by if/ over the complementary of a certain compact interval. We follow 
[121 and write if/ as if/ - if/ + C, where if/ e '^^(U, W), and verifies if/ - -C over for N large enough. Using the 
technique developed in (12| based on (7) and Poincare inequality, we have 

Var [Tri/A(XivZ^)] = Var [Jr iff a n:L*^]] = & (-i^] , 

E[Jrif/{-LMl.*f^)]^M f v>(A)d^iN(A) + (^| 
AsE[Tri//(ZjvZp] = CM+ E [Tri/i'(ZjvZ;;r5]' *is leads to E [Tn//(ZAfZp] = ^ (^). As 

E[(Tri/A(ZjvZ;v))^l = (E[Tn/'(ZjvZ;,)])^ + Var[Tn/'(ZjvZ^)] (30) 

we finally obtain that (29) holds for / = 1 . 

We now assume that I l29t holds until the order / — 1 for each function of S^'^ (R, [R) vanishing on for N large 
enough and constant over the complementary of a compact interval. We consider such a function if/ and evaluate 
the behaviour of the 2/-th order moment of Tri//^(ZjvZ^). We have 

E[(Tn/.(ZivZ;;,))''] = (E[(Tn//(ZwZ;v))'])'+Var[(Tn/.(ZivZ^))'] . (31) 

The first term of the r.h.s of (31) can be upperbounded as follows 

(e [(Tn/^(ZwZ^))'])' < E [(Trt//(ZA,Z^))'] E [(Tri/.(Z^Z;;,))'"-"] = (^] , 

using that (29) holds until the order / - 1. The second term of the righthandside of (3D can be evaluated using the 
Poincare inequality. Using that the partial derivative of Tn/'(ZAfZ^) w.r.t. VKj.j.jv and Wi,j,jv are equal respectively 
to eJZ^i/''(ZjvZ^)e; and ej y/' {l,i^l,*^]'I,i^e j , we immediately obtain that 











= -MC + 


(^) 









Var[(Tn//(ZjvZ;v))'l ^ CE 



-^Tr (i/a'(ZjvZ^)2ZwZ^) (Tn/'(ZjvZ^))^' ^ 
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Using Holder's inequality, we get immediately that 



Since the function A yj'{X)^X belongs to 
holds that 





/ 


<c 


E 







(E[(Trv^{Z^Z;;,))^']) ' 



(32) 



and has a support disjoint from for N large enough, it 



^Tr([i//'(2wZ;,)]'2iv2^) '< We ^-Tr ([v^'(Z^,Z^)]'Xjv2:;v) T^Tr ([v^'ffivZpl'XjvZ;;,) 



1 

N 
1 

Plugging the previous estimates into (32), we get 



2(/-l) 



Var[(Trv.(ZwZ;;,))'] < ^ (e [(Trt/.(Zjv2^)f ]) ■ 
Define xjv = E |(Tri//(Z]vZ^))^'j and un - iV^'xjv- From (31), we have the inequalities xn < -p-J^jy' + and 



i-i 
I 



um £ Ci + C2. We claim that the sequence [un) is bounded. If this is not the case, it exists a subsequence u^j^ 
extracted from un which converges torwards +00. However, the inequality + > 1 must holds for AT large 

enough. As w^j^ — +cx>, this leads to a contradiction. Therefore, wjv is bounded and xjv £ for iV large enough 
This proves Lemma|3l 



□ 



3.2 Evaluation of the escape probability of the id) yt,jv) k=i,...,N 
In this section, we will prove the following result. 

Proposition 3. Assume assumptions A- 1 toA-6 hold. For each I e H, it holds that 



We foUow the same approach than in Section lSJi and first prove that the {(hk,N) k=i,...,M satisfy a property similar 
to (7). For this, we study the behaviour of the Stieltjes transform /2]v(z) of the distribution -pZ^i ^A/tjv defined by 

1 ,-1 
"jv(.z) = — Tr^fliv-zIJ . 

and use Lemma [T] as well as the inverse Stieltjes transform formula (3). Our starting point is the following result 
showing that the empirical eigenvalue distribution of fljv is very similar to the distribution of the eigenvalues of 
ZjvZ^. The following auxiliary result will be useful. 

Lemma 4. Assume assumptions A- 1 to A-6 hold. It holds that 

^[m'^{z)]-m'f^[z)^^^ (33) 



where tjv is analytic on C\IR and can be upperbounded byPi (|z|)P2 1 |jm(z)| ) ""^ " 
Proof: The proof is given in AppendLx l5.II □ 
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We now prove the fundamental following result. 
Lemma 5. Assume assumptions A- 1 toA-6hold. For each ze 



E[«iv(z)]= + — — + — rr?-. 



withKM afinite signed measure carried by S^M such thatXM^lx^ N'-^-q n^] ~ ^fa^l - ^•■■■< Qn? andr^ aholomorphic 
function on C\R satisfying 



|rjv(z)|<Pi(|z|)P2 



1 



JIm(z)h 

with Pi , P2 two polynomials with positive coefficients independent ofN. 

Proof: Using that Iliv is a rank 1 perturbation of Aw, we obtain immediately that 



1 cr^CNin'^iz] 



M 1 + a^CNrfiMiz) 
Therefore, for z e C\U, it holds that 

a^CNtn'^iz) 



E[n]v(z)] = E[mivU)]--J-E 
M 



l + cr^cjvmjvCz) 



(34) 



We first establish that 

a^CNm'j^(z) 



o'^c^m'iz) cjv 

~ +^r„(z), (35) 



l + cr2civmjv(z) N 



l + a^c^misiiz) 

where rN^z) is holomorphic on C\IR and satisfies |r]v(z)| < Pi(|z|)P2 For this, we write 

a^CNm'^iz) a^CNm'j^iz) 
l + a^CNtfiNiz) 1 + cr^cjvm]v(z) 

a^CNim'j^iz) - m^(z)) (cr^cjv)^ (mjv(z)(m^(z) - m^(z)) + m^(z)(mjv(z) - mjv(z))) 

il + a^CNrhN{z)){l + cr^CNmNiz)) {1 + a^CNmN{z]){l + a^CNmf^[{z]) 



(36) 



In order to study the expectation of this expression, we use i fTTi and (15). Moreover, (6) and a straightforward 
application of the Poincare inequality to mjv(z) considered for z fixed as a function of the entries of leads 
immediately to 



E|mjv(z)-mjv(z)|2< — :rPi (|z|)P2 



1 I 1 



Af2 i|Im(z)h 

for some polynomials Pi , P2 with positive coefficients and independent of N. Therefore, 

E|mjv(z)-mjv(z)|< ^[Pi(|z|) + P2 ^ 



N[ Ulm(z)|; 
Applying also Poincare inequality to bound Var[m^(z)], together with LemmaH) we get 



E|m;,(z)-m;,(z)p<^Pi(|z|)P2f ^ 



Therefore, it holds that 



lllm(z)| 



Elm' (z)-m' (z)|< ^ Pi(|z|) + P2 ^ 



N\ Ulm(z)|, 

Using I mjv (z) | < I , | m'{z) | < ^ 2 - as well as (11) and (15) , we eventually get from (36) that 



E 



a'^CNtn'j^iz) a^CNm'^iz) 



1 + a'^CNifiNiz) l + a-^c^^mpjiz) 



1 I 1 

<-Pl (|Z|)P2 



N Ulm(z)| 
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This immediately implies (35). Now define the function h^iz) by 



1 + a^CfftriNiz) ' 



This function coincides with the Stieltjes transform of a signed measure ?Cjv satisfying the conditions of Lemma 

5- where is a compact containing 



\5\ Using QOJ, we obtain that \hM{z)\ < 2a^CN\m'^(z)\. As \m'^{z)\ < ^(^^^j; 



5^jv, it holds that |^jv(z)| < C- 



Tj. Using Theorem 4.3 in O, we obtain that hMiz) is the Stieltjes transform of a 



finite signed measure Kiv, the support of which is the set of singular points of h^iz), i.e. S^n. In order to evaluate 
KNi[x~ j^,x'^ j^]), we use the inverse Stieltjes transform formula, 



(/ /ijv(x+/y)dx . 



It is clear that /ijv(x + iy) - ^'°gti+°' ovmArCx+iy)) ^ ^j^gj-g ^.j^g complex logarithm corresponds to the principal deter- 
mination defined on C\[R~. We note that i fTO) justifies the use of the principal determination. Therefore, 

/ ^ /ijv(x+jy)dx = log|l + o-^Cjvmjv(x^^^+iy)j-log|l + CT^Cjvm]vU;^^jv+'y)]- 



When y — 0, this converges towards log(l + (7^CArmAr(x^j^)) -log(l + ct^ cat mAf(x^ j^)), a real quantity because x' 
and x"^ ^ belong to dS^N- This shows that jcat ([x~ ^, x"^ ^]) = 0. Consequently, 



a^CNm'j^(z) 
1 + a^CNrhNiz] 



-L 



dKNW CNVMiz) 

^„ A-z ^ N ' 



where rNiz) is holomorphic on C\IR such that |rjv(z)| < Pid^DPaljj^^j)]-)- Lemma [5] follows immediately from OH . 

□ 



We now handle the proof of Proposition|3l Although certain steps of the present proof are similar to the proof 
of Proposition|2] more work is needed because matrix {In considered as a function of the entries of Wjv is more 
^. We still consider function y/Q e (U, defined by (28) and remark that 



complicated than Z]vZ 



P(^2,jv) < P(Tri//o(«w) > 1) < E [(Tiy/oiClN)f 



for / e Kl. In order to establish Proposition^ it is therefore sufficient to prove thatE[(Tn/'o(«N))^'] for 
each integer /. For this, we still use the Poincare inequality. However, in contrast with the context of Proposition^ 
the entries of iljv, considered as functions of the real and imaginary parts of the entries of Wjv , are not continuously 
differentiable on IR^*^^ because function Wjv ^k,N is not differentiable at points for which eigenvalue A^ jv is 
multiple. The use of Poincare inequality has therefore to be justified carefully. The following useful lemma is 
proved in the appendix. 



Lemma 6. Assume assumptions A- 1 to A-6 hold. Let ifr be a function of '^^{U.,U). Then, Tiij/{ClN), considered 
as a function of the real and imaginary parts of the entries ofW^, is continuously differentiable. Moreover, if the 
have multiplicity 1, it holds that 



eigenvalues ofL^^ 



,—T:r(if'[ClN)]\^ — 



M 

1=1 



l,N 



(37) 



where ft/^jv represents the orthogonal projection matrix on the 1-dimensionnal eigenspace associated to the eigen- 
value X in o/ZjvZj^. 



We will also need that 
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Lemma 7. For each integer p>0, it holds that 



supE[||WwW;,|| 

N 



< +00, 



a property also established in the appendix. We now prove the following result. 

Lemma 8. Assume assumptions A- 1 toA-6hold. For all function y/ e {U,U) constant over the complementary of 
a compact interval and which vanishes on the support Sfi^q of p n for all N large enough, it holds that 



(38) 



for each / e N. 

Proof: As previously, we prove Lemma[8]by induction on I. We first consider the case I - I, and consider a function 
Iff as above, and denote by C the constant value taken by y/ over the complementary of a certain compact interval, 
and by y/ the function of '^^{U.,U) defined by y/iA) - y/{X} - C, which, of course, is equal to -C on 5^n- Using 
Lemma[T]and LemmaO we obtain 



E 



^Trv>(njv)] = yrWdpNW + ^j^ V>(A)dKjv(A) + ^[-^ 



Using that ^nHx' j^tX'^ ^]) - for each q - 1, . . . , Qjv, we get that f^^y/{X)dKisi{X) - and that 



(39) 



E 



^Jiy/iClN) 
M 



Therefore, it holds that 



E 



-!-Tn//(njv) 
M 



Moreover, we prove the following lemma. 

Lemma 9. Assume assumptions A- 1 to A-6 hold. It holds that 



1 



Var 



^Tr(v/(niv)) 



(40) 



Proof: We first note that, considered as a function of (Re(W;j-,jv),Irn(Wi,j,jv))i<i<M,i<j<iv. function j^Try/{ClN) is 
continuously differentiable by LemmalU Therefore, function j^Jj:y/{ClN) is continuously differentiable as well. It 
is thus possible to use the Poincare inequality to evaluate the lefthandside of HOl . Furthermore, as the probability 
that the eigenvalues (Aj; ]v) k=i,...,M have multiplicity one is equal to 1, it is sufficient to evaluate the partial deriva- 
tives of function j^Jiy/iflN] when is such that the (A,jt,iv)fc=i,...,M have multiplicity 1. As the derivative of y/ 



coincides with yr , 137) and Poincare inequality lead to 



Var 



-Tr(i//(fiw)) 



C 

< ^E 



1 

— Tr 

M 



M 



ZwZ^^|[i/A'(njv)],,/| nz,iv 
1=1 



or equivalently. 



We claim that 



Var 



-Ti[yf{nN)) 



c 



M 



^ /—I 



|[/(njv)]/,,| <{iyf'{nN)f)ii. 

Indeed, if (Vfc,iv)fc=i,...,M represent the eigenvectors of Q, then 



(41) 



M 



[i//'(niv)]„= Y.V'i'^k,N)\eJvk,N\^ 

k=l 
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As V|fc,ivl - 1, Jensen's inequality yields to (41). Therefore, it holds that 



Var 



-Tr(i/.(n^v)) 



Af2 



I M 



As supjY IIBjvB^li < +00, we get using lemma[7]that 



supE[||ZivZ;^||P] <+oo. 

JV 

We remark that || 2jvX^|| < t2 +eon the set j^,and write the righthandside of 02) as 

M 



c ^ 



It holds that 



M 



M: 



1=1 



M 



Tr(i//'(njv))^ 



Function tfr'^ belongs to '^^{M.,U) and vanishes on^^v- Therefore, lemma[5]implies that E [-pTr (i//'(n]v)) 
(see Eq. (39)). Moreover, as -pTr [y/'iQ.N)^) < sup;^ i//'(A)^ < C, we have 



Y M 



which is itself upperbounded by 



for each integer p. This completes the proof of lemma[9l 

Assume that (38) holds until integer / - 1. We write as previously that 

E [ {TTy/{nN)f ] = (e [ (Tr v(njv)) ' ] )^ + Var [ [Try/iClN))' 
The Cauchy-Schwarz inequality leads immediately to 

(e [(Tri/^Cliiv))'])' < E [(Tri//(niv))'] E [(Tn//(niv))""'] = ^("j^) ■ 
As for the second term of the r.h.s. of (43) . we use Poincare inequality and Holder's inequality to obtain 

Var[(Tri//(niv))'] <CE 



?r;-n 1 *f . 9 
(TrV'(nw))'""'^ — ^ Xk.N{[V'lilN)h.k) 
™ k=l 



< c 



1 ^ - . 2 

™ k=l 



(E[(Trt//(n^))'']) ' 



Jensen's inequality leads again to 



M 



1 - ^ 2 

■TZ E ^k,N{[V'i^N)]k.k) 
™ fc=l 



< E 



ZNZ*j^\\—Jr:y/'[nNy 



We write again that 



l.N^*N\\—Jry/'(£lNr 



M 



IZiv 2:^11— Tn//'(njv)2ll<fi,, 



M 



+ E 



ZatZ^II— Tri/A'(niv)2 



M 
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and obtain as previously that 



Xiv2:«l|— Tri//'(njv)' 



M 



<C 



—Trnf'inMf 



1/2 



But, applying Cauchy-Schwarz inequality as in (44) to E|Tri//'(njv)^|' leads to E| j^Tri//^'(fi]v)^|' -<^[-^]- Gather- 
ing all the previous inequalities, we find that 

E [[TryfinM)f] < ^ (e [ (Trv^(n^))'']) ^^[w]' 

and in the same way as in the proof of Proposition |2] we obtain E |(Tn/'(njv))^'] = (i^j- "^^^^ concludes the 
proof of Lemma|8l □ 



3.3 End of the proof of theoremU] 

We now complete the proof of Theorem|5]when function 9 — ► a(0) is given by 



for Be [-jT,7i]. We recall that Sff is defined by 



he" 



j(M-\)e 



where (ifi,iv)i=i,2 are defined by l|24) and | |25) , and that tgc^ < xn where xn - det (/)(Zjv2^)det (p{Clisi). We first give 
a useful lemma which appears as a straighforward consequence of the evaluations of Section l23l 

LemmalO. Assume assumptions A- 1 toA-6hold. ForeachN, itholdsthat 

sup ||Tiv(z)||<C, 



sup 



1 



sup 



\ + o'^CNmNiz) 



<C, 



and for N large enough, we have 



sup ||Qjv(z)||xjvS C, 



sup 



Xn 



sup 



1 + a'^CNmj^(z) 



< C, 



1 + (J^Cf^mi^{z) 



XN<C. 



We consider the set 



, 2{k-Dn , n 



and remark that for each 6 e {-n,n] and for each N, there exists 8^ e such that \6 - Sjvl S For each e 
[-H', :/rl , it holds that 



f]N{0) - r]N{0) = [flNiB] - flNiBN)] + [flNiBN) - r]N{BN)] + [riNiBN) - 7?jv(0)] . 



(45) 
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It is easy to check that the third term of the r.h.s. of (45) satisfies 



sup |77jv(0jv)-7?iv(0)|<2 sup ||a(0)-a(0jv)|| 



0e [-71,71] 



6e[-7r,7r] 



(46) 



In order to evaluate the behaviour of the supremum over 9 of the first term of the r.h.s. of 145) , we prove that for 
each a > 0, 

/ 1 



sup |TjAf(0)-fjjv(0iv)l > a 
^6E[-7t,7r| ) 

where /3 > 0. We first remark that for each / £ N, it holds that 



^0£|-7I,7t] 







1 


1 > a 


< P 


sup 1 






^eE[-7t,7r] 



sup \f]N(d)-f]N(dN)\\-<^> a 



1 

<-E 
a' 



sup |f?iv(0)-j?jv(0]v)| 

0£[-7r,7r] 



+ P(<?JV) 



1 



Moreover, 



|?7jv(0)-j?jv(0]v)|'%,<cf 



(a(0)-a(0jv))*Qjv(z) 



w'{z) 



Lemma[TO]and the inequality l^c ^ Xn imply that 



sup ||Q]v(z)|| 



1 + a^Cf^rhf^iz) 



■\gc < C 



a(0]v) 



Idzl. 



for some constant term C. Inequality (46) thus implies that 



sup 

8£[-n,n] 



(a(0)-a(0iv))*QwU) 



1 + a^CNrriNiz) 



thus showing that 



sup |7jjv(0)-T)jv(0Ar)| > a 

^0£[-7t,7r] 

for each integer I. Borel-Cantelli's lemma eventually implies that 



a(0iv) 



1 



(47) 



sup |77jv(0)-??jv(0jv)| ^0 
ee[-7i,7r] 

almost surely. 

We finally study the supremum of the second term of (45) . We denote by Vk,N the elements of d^. Let a > 0, 
then 

\ I 

sup |f7jv(0iv)-T7iv(0iv)| > a £ P sup |f7jv(Vfc,iv) - J7jv(vi;,iv)| > a 

\ee[-7r,7i] / \fc=l,...,iV2 y 



^ P (|77Af(Vfc,jv) - ?7iv(Vfc,jv) I > a) 
fc=i 

E [P({|'7w(Vfc,iv)-J7iv(vfc,jv)|>a;}n^w)]+^ 
fc=i 



1 



for each integer /. We now introduce in the above term the regularization term xn - det (p{l.isil.*^) det (p{flN) 
defined in (27) . As jjv is equal to 1 on <f^, it holds that 

P {{\flNiVk,N) - VN(Vk.N) I > a} n .g"^) = P ({|7?]v(Vfc,jv) - VNiVk,N) \x%>a}n S'^) 

< P{\flN{Vk,N) - '7iv(Vfc,Jv) I > 



y2l 



E I ('7Jv(Vfc,iv) - r]N (Vfc,jv)) 



,21 
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The introduction of jjv is in part motivated by the observation that the moments of 77]v(vj;,iv)X^- are finite. More- 
over, it holds that 



E I (jyjV (Vfc.Af) - ??Af (Vfc.iv)) X jv 



2 |2/ 



a(Vfc,iv)* QivU) 



-Tiv(z) 



1 + a^cisiniMiz) 



a(vj;,jv)X 



N 



21 



\dz\ 



(48) 



In order to complete the proof of Theorem|5l we establish the following proposition. 



Proposition 4. Assume assumptions A- 1 to A-6 hold. TfCajvJweM is sequence of deterministic vectors satisfying 
liaivli - 1, then, for each integer!, it holds that 



sup E 



■Tiv(z)- 





21- 


jajvlw 





c 

N' 



1 + a^CNrhisiiz) 1 + a^CNmN^z) i 

where the constant C does not depend on the sequence (ajv) . 

Proof: In order to shorten the notations, we denote by gwCz) and gwCz) the functions defined by 

w'j^iz) 



(49) 



and 



gNiz] = a^Qjv(z)aiv 



gN(z] = a^Tjv(z)aA, 



1 + a^CNniNiz] ' 



w'j^iz) 
l + cr2cjvm]v(z) 



In order to evaluate E|gjv(z) - gNiz)x%f'' , we use the Poincare inequality. For this, we first state the following 
lemma proved in the appendix. We recall that if H a hermitian matrix with a spectral decomposition H - Z/ 7/X/X*, 
its adjoint (i.e. the transpose of its cofactor matrix) denoted by adj(H) is given by adj(H) - Zz(nfc7^z7fc)x;Xp When 
H is invertible, adj(H) - det (H)H~^. Next, we state the following lemma proved in the appendix. 

Lemma 11. Assume assumptions A- 1 to A-6 hold. Considered as functions of the real and imaginary parts of the 
entries ofWN, functions det (/)(ZivZ^) and det (p{ClN) belongto'^^iM^'^^), and their partial derivatives w.r.t. VKi,j,iv 
denoted by 



[Di];j,jv:= 
[D2]/,j,]v:= 



d 



dWi 



>.j,N 



{det(/)(ZivZ^)}' 
{det(t>iQ.N)}, 



are given almost surely by 



[Dilij.jv = e*Z^adj(^(Z]vZ^))(/)'(ZivZ^)e,-, 



VD2\i.j,N- 



M 



Z [adj ((/)(niv)) </)'(niv)] „ ft/, 



N 



(50) 
(51) 



J' 



If we denote by s^i^n cind si2,N the events defined by 



sdi^N = {3fc : Xk,N ^^e}r\ {Xl,N, ^M.N £ supp((^)} , 

si2,N = {3fc : aicN^Sr^.] n {wi.iv, ■ ■ .,&m,n e supp((/))} . 
then [Dilj ^- jv - onsi^j^ and [Ti2]i,j,N - onsi^j^. 

We now establish l l49t by induction on /, and first consider the case / = 1. We write the second moment of 

igNiz] - gN[z])x]^ as 



E I (giv(z) - gjv(z))x^|^ = |E ((g]v(z) - giv(z))x^) 1% Var ((gjv(z) - giv(z))x^) . 



2 1|2 
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We evaluate Var [ (gAf (z) - gN{z))x%] using the Poincare inequality and get 



Var[(giv(z)-gwU))X^]<— 

i.j 



It is clear that 



We verify that 



so that 



i.J.N 



dgNiz) 



W 



i.j.N 



+ 2E 



\gNiz) - gNiz)\ 



i.j.N 



dWij^N dWij^N l + cr^cjvmivU) 



+ a^Q]v{z)ajv- 



^^iz) \ 



dWij^N 1 1 + a^CNniNiz) j ' 



dQNiz) , , 

a„rrTTT ajv - -a^QNiz)ei eyZ^Qjv(z)ajv, 



E 



dWij,N 

2 



dWi 



i.j.N 



aMQiv(z)Qiv(.z)*ajvatQjv(z)Zjv2tQiv(z)*aiv 



As Xn ^ implies that Am,n - ll^^jvZ^H < + ^e, LemmafTOlimplies that 

sup ;(;^a]^Q]v(z)QAr(z)*aiva^Q]v(z)Z]v2]^QivU)*a]v< C. 



(52) 



Using again Lemma [TO] we get that 



sup Xn 



Wj^iz) 



1 + a^CN>nNiz) 



E 



dQNiz) 

'dWi 



-ajv 



i.j.N 



< c. 



We obtain similarly that 



sup I^|a^Qjv(z)ajvl^^ 



5W,-,j,jv 1 1 + a^CNmNiz] ) 



C 



(53) 



(54) 



The same conclusions hold when the derivatives w.r.t. variables VVj.j.jv are considered. This shows that the first 
term of the r.h.s. of (52) isa& [-^] term. We now evaluate the behaviour of the second term of the r.h.s. of (52) , and 
establish that 



sup E 



zei3S?„ 



for each integer p. We express g^J^ as 2xNg^^- Lemma [10] implies that sup^^g^^ x^|gjv(z) - gjv(z)|2 < C. 
Therefore, it is sufficient to check that 



\gNiz)-gNiz)\^Y. 

'.J 





2i 




dWi^j^N 




\NP] 



(55) 



dXN 



Wi 



i,j,N 



1 

JVP 



for each integer p. ^^'^^ can be written as 



dx 



N 



dWi 



It holds that 

E 



i.j.N 



Y^\lDi]ij,Ndet(p{ClN))\ 

'.j 



[Di ] i,;,jvdet (P iClN) + [T>2 ] i,j,ivdet ^(ZivZ^) ■ 



= E [det </)(nw)^Tr (Xjv2:w'^'(2:iv2:w)^adj ((/)(2ivZ]v*))^) M.i 



Moreover, we can write 

Tr [zN-L*j^cl)'i-LN:^N^^adi ((^.(XjvXjv*))^) ^Y.hN(p'ihNf fl (pihNf 



;Tr(ZN2:>'(ZwZ;v)^). 
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because (p{A) < 1 on K. Therefore, it holds that 



E >:i|D,l,:.;jvdet(/)(nA/)r <E[detil)[nNfTi[i:Nl.lr(p'i^N^*Nf]M,J^CNP[s^i,N), 
because det cp{ClN)) < 1, and Tr (ZjvZ^(/)'(Zjv2;^)^) < CN on s^i^n- As ^i,iv <§'i,jv, proposition|2]implies that 

E 



X|[Di]ij,ivdet0(niv))| 



1 



for each integer p. Using similar calculations and Proposition|3] we obtain that 



E|[D2],-,;,wdet(/)(ZjvZ;v))| 



1 

'np 



for each integer p. This completes the proof of (55\ and establishes that 



sup Var [(gjv(z) - g]v(z))i^) ^^{j;^ 



In order to evaluate the term |E [(gjv(z) - gN{z))x%] f, we also need the following auxilliary lemma proved in the 
appendix. 



Lemma 12. Assume assumptions A- 1 toA-6hold. It holds that 

sup |E [a^QAr(z)ajvXiv - a^Tiv(z)ajv] | = G 



I 1 



1 



sup |E[miv(z)xjv-miv(z)]| 

sup |E[m;vt^)iJv-m;v(^)]h^h^ ■ 

We express (gjv(z:) - giv(z))x^ as ;6i,iv(z) + Pz.nIz) where 

/5i,iv(z) = Ijv (a^QAf(z)aiv-a^TAf(z)aivJ 5 

Pz.Niz] = x^a;yT]v(z)aivl — 

and establish that 



(56) 
(57) 
(58) 



and 



/32,iv(z) = x2^a^Tiv(z)aiv ,^2^ r , " T" ■ 

\. 1 + CT'^C]vmAf(z) l + CT'^Cjv?Wjv(z)j 



sup E\Pl,Nf = < 



and sup \E\p2,N\ 



ZEdSly 



1 



(59) 



Using Lemma [TOl (59) for /3i,jv will be established if we show that 

sup E|xjv(a^QAr(z)ajv-a^Tjv(z)aAf)| 



For this, we write that 



^\XNia*NQNiz]aN - a^Tjv(z)aAf) 1^ = Var(Xiva^Qjv(z)aAf) + |E (xAr(a^QA,(z)ajv - a;;rTiv(z)ajv)) f . 

The above calculations prove that sup^^^^^ Var[X]va^QAr(z)aAr] - 6 (-^j, while (56) and 1 - E(xjv) = for each 

p imply that 



E [ X jv (a^Qjv (z) aw - a^Tw (z) ajv) ] 



>fJ_ 
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This completes the proof of (59) for Pi,n- In order to show (59) for Pz.n, we first remark that by Lemma [TOl 
|a^T]v(z)ajvl is uniformly bounded on dS^y, and write that 



m'{z) 



m'j^iz) 



1 + a^CNifiNiz) 1 + (T^cjvmjvCz) j 



1. 



or equivalently that 



w'j^iz] 



w'j^iz) 



l + CT^cjvmjvtz) 1 + a^CMiriMiz) j 
o'^cjvX^ irfiNiz) - niNiz)] + 2z(t^cnxIi ('^iv^^^ ~ '"w^^)) 



■iCT^CNfa^a-CN) 



- cr^ cjv 0-^(1 -cjv)- 



Xn 



(1 + a^CNmN{z)){l + a^CNniMiz)) 



[mAr(z)xjv(m^(z) - m^(z)) - m^(z)xAr(mAf(z) - mjv(z))] 



(1 + cr2c^^^(2))(i + a'^cj^mj^iz)) 
The Poincare inequality and Lemma[T2]imply that 



Xjv(m^(z)-m^(z))]. 



sup E|xjv(mjv(z)- mjv(z))|^ = 



and 



Eq. (59) follows immediately from 



sup E|xjv(m^(z)-m^(z))| 

z£i3S?„ 



1 

iV2 



sup 



liV 



(1 + a^CNtfiMiz)) (1 + a^CNrriMiz)) 



for some deterministic constant C (see Lemma[TO). This completes the proof of (49) for / = 1. 
We now assume that (49) holds until integer / - 1 and write that 

^\X% [gNiz) - gNiz)]f' = |e [(i^(giv(z) - giv(z)))'] I' + Var [(x^(gw(z) - gjv(z)))' 
The Cauchy-Schwarz inequality implies that 

|e [ [XNigNiz] - giv(z)))'] I' < E |x^(giv(z) - giv(z)) f E |x^(giv(z) - giv(z)) | 



,2{l-l) 



and shows that 

The Poincare inequality gives 

Var[(i2^(gjv(z)-giv(z)))' 

Finally, (53) and (54) imply that 



sup E(;i;^(giv(z)-giv(z))) =( 



1 

AT' 



N 



XN^gNiz)- gNiz)) 

gjv(z)-giv(z)) 



2(;-l) -,4 V 



gNiz) 



+ ^E 



N 



sup iwL 



Z£dSiy 



I. J 



dgNiz) 



dWi 



i.j.N 



dgN(z) 



|2' ^4Z-2y- 



<C, 



dWij,K I 

'^;t:]v 



gNiz) 



for some deterministic constant C. Therefore, the supremum over z £ of first term of the r.h.s. of (60) is a 

^ [■j^j- Moreover, it can be shown as in the case I - 1 that the supremum over z £ dS^y of the second term of the 

righthandside of (60) is a ^ iw) each integer p. This completes the proof of Proposition|4]and of the uniform 
consistency of estimator jyjv (0) . □ 
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4 Consistency of the angular estimates 



We now adress the consistency of the DoA estimates defined as the local minima of function 6 fjNid)- For 
this, we assume that the number of sources K is fixed, i.e. that K does not scale with N. In other words, model 
l.N -Bn + Wjv corresponds to a finite rank perturbation of the complex Gaussian i.i.d. matrix W^f . 

Remark 2. In this context, it is possible to derive in a simpler way than above an alternative consistent estimator, 
say 9 fjN.spikeiS) of function 6 ^ rjNid). This estimator is obtained by assuming from the very beginning thatK 
is fixed, and is based on the recent work ofBenaych & Nadakuditi 13]; see 122] for more details. However, as shown 
in \22} , estimator f]!^ [6) always leads in practice to the same performance as fiN^spikeid) ifjj « 1 (typical value 
in !22l), but outperforms fjNspij^eiS) for greater values of jj (typical value ^ in f22ll ). Therefore, the use of estimator 
f] N [8) appears in practice more relevant than f]f4 spike W ■ 

In order to define the estimators of 0i , . . . , 0^ properly, we consider K disjoint intervals Ii,...,If^, such that 
Qfc e Ik< and define for each k the estimator Qfc,N of (^k by Qk,N - argminge/j. \ fi]^{d)\. We prove the following result. 

Proposition 5. For k-\,...,K, with probability one, 

mBk.N-Ok) — -0. 

Af— oo 

In order to establish the proposition, we follow a classical approach initiated by Hannan [15] to study sinusoid 
frequency estimates. For this, we first recall the following useful lemma. 

Lemma 13. Let [um) a real-valued sequence of a compact subset o/ (-0.5, 0.5], and converging to a as M ^ oo. 
Define quia m) - TjZtli e~'2'i'=«M r^Q. ^ g or if a - QandM\aM\ oo, thenquiocM) 0. If a - andMaM '■ 

peU, then quiocM) ^ ■ 

We denote by A the matrix A(0) corresponding the true angles = {8\,...,QkV . It is clear that rj^iO) = 1 - 
a(0)*A(A*A)"iA*a(0). By the very definition of Bjc,n, \f]N{dk,N)\ ^ |jjiv(0fc)l- From Q) and the equality r]N{0k) = 0, 
we have |^jv(0fc,iv)l ^ w.p.l., as JV ^ oo. Consequently, 

\VNidk,N)\ S \riNidk,N)-flN{dk,N)\ + \flN{dk,N)\ 

< sup |7?jv(0)-7?Jv(0)I + I??]v(6'a:,]v)I 
0£[-7r,jr] 

0. (60) 

From Lemma [13] (A*A)~^ converges to Ijf as AT ^ oo. Since (0jt,iv) is bounded, we can extract a converging subse- 
quence {Ok,(p(N))- Lef oiN - dk.cpiN) ~ dk- From Lemma[T3] if ajv ^ as Af ^ oo, then 

a(0fc,^(iV))*A(A*A)-iA*a(6ifc,^(jv)) 0, (61) 

and thus r]Nidk,ip{N)) ^ 1. a contradiction with | |60] . This implies that the whole sequence (0fc,iv) converges tor- 
wards 0jt. If ^l6*fc,Af~'5fcl is not bounded, we can extract a subsequence such that N|0fc,0(jv) -0fcl ^ +oo and Lemma 
[T3] again implies that l l6Ti holds, a contradiction. N\6]c,n ~ Sk\ is thus bounded, and we consider a subsequence 
such that Ni6]c,ip{N) - Ok) ^ P where pe [-n,n]. From Lemma [T3] if /3 0, we get 

flwiN) i9k,(p(N)) ' 1 ~ 

N^oo 

which is again in contradiction with (60). Therefore, (} - and aU converging subsequences of (A/^|0fc,,/,(Ar) -dk\) 
converge to 0, which of course implies that the whole sequence (A^|0fc,Af ~ ^ifcl) converges to 0. We finally end up 
with NidicN -6ic)^0 w.p.l., as A/^ ^ oo. 
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5 Appendix 

5. 1 Proof of Lemmall} estimate of E [ m'^ (z) ] 

We first give the following useful technical result. Its proof, based on Poincare's inequality, is elementary and 
therefore omitted. 

Lemma 14. Let (Mjv(z)) a sequence of deterministic complex M x M matrix-valued functions defined on C\U such 
that 

||MAr(z)|i<Pi(|z|)P2(|Im(z)r^). 

Then, 

Var 
Var 



-^TrQ]v(z)Mjv(z) 



<^Pi(|z|)P2(|Im(z)r^), 



N2 



■^TrZ^Qiv(z)Mjv(z) 



<-^Pi(|z|)P2(|Im(z)r^). 



Moreover, the same results still hold when Q]v(z) is replaced foj/Qjv(z)^. 
We are now in position to establish Lemma|4l We have to establish that 

|E [m^(z)] - m'j^(z)\ < -^Pi(|z|)P2 (|Im(z)r^). (62) 

For clarity, we recall results from fW], f23l and fT4l, on which the proof heavily relies. We have first to introduce 
some new notations extensively used in |10|, |23| and 1 14|. We define 5jv(z) - acf^mj^iz) - cr-^Tr (Tjv(z)), as well 
as 5jv(z) - Sjsiiz) - "^^'J^'^'' which coincides with the Stieltjes transform of finite measures cjvjUn + (1 - cjv)5o. In the 
following, matrix T]v(z) is defined by 

tjv(z)= -z(1 + (T(5jv(z))Im + 



l + cr5jv(z) 

and is related to 5jv(z) through the equation Sn(z] - o--^Tr(f jv(z)) (cf [13], [23]). We also remark that matrix Tjv(z) 
can be written as 

Tiv(z) = -z(l + (t5jv(z))Im + ^ , 

[ l + cr(5iv(z)J 

and that wn{z) coincides with z(l + (t5jv(z))(1 + (t5jv(z)). We also denote QAf(z) the resolvent of matrix Z^Zw, i.e. 

Qjv(z) = (Z^Zjv - zljv)"^ 
and define a^iz] - E [■^TrQjv(z)], ajv(z) = E [-^TrQArlz)], and the matrices 



Rjv(z) = |^-z(l + CTajv(z))lM + 
Riv(z) = -z(l + (Tajv(z))Ijv + 



1 + aaN(z) 

B^Bjv 
1 + aa^iz) 



It is shown in [lOJ and |23| that the entries of QAr(z) (resp. QAf(z)) have the same behaviour as the entries of 
Rjv(z) and Tjv(z) (resp. of Riv(z) and TAr(z)). It is also useful to recall that \aN{z)\, |aiv(z)|, \-z{l + craNiz))\~^, 
|-z(l + (Tajv(z))r\ ||Tjv(z)||, ||Tjv(z)||, ||Riv(z)|l and ||RAr(z)|| are bounded on C\IR by Pi(|z|)P2(|Im(z)ri). We re- 
mark that our new notations are symetrical w.r.t. the substitution Zjv ^ Z]^, and are easier to use in the forthcom- 
ing calculations. 

We first notice that l l62l is equivalent to 

\a'^{z)-S'^{z)\ < -^Pi(|z|)P2 (llm(z)r^) - (63) 
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In order to prove (63), we first show that 

\a't^iz)- ^TrR'^{z)\< ^Pi{\z\)P2{m[z)\-^) (64) 

and deduce from this that (63) holds. Using results on the behaviour of ajv(z) - -^TrRjvCz) established in [10], [23l 
and [M], we first establish that l|64) holds. For this, we recall the following lemma. 

Lemma 15 ([10], |23l proof of Prop.6]). Forz e C\R, it holds that 

E [Qjv(z)] = RjvU) + Aw(z)Rjv(z) + [^Tr Ajv(z)] E [Qjv(z)] Riv(z) (65) 
where An{z) is given byA^iz) - Ai,jv('Z) + A2,Af(z) + A3,jv(z) with 

Ai,iv(z) = E [Qjv(z)ZwZ^-^Tr (Qjv(z) - E [Qiv(z)])l , 

1 + (7q:jv(z) ^ Af J 

^2 



A2,iv(z) = E [(Qiv(z) - E [Qjv(z)]) ■^Tr2;,QAr(z)Bjvl , 

l + (TaA/(z) ^ _/V J 



A3,iv(z) = - , ^^. E [Qjv(z)] E [ -^Tr (Qiv(z) - E [Qjv(z)]) ■^TrZ;^Qjv(z)Bjvl . 

7fMjv(z) «s a sequence of deterministic complex matrix-valued functions defined on CW such that 

||Miv(z)|l<Pi(|z|)P2(|Im(z)ri), 
then j^Tx: [\i^i^{z)Mp]{z)) , for i - 1,2,3, is bounded by ^Vi{\z\)¥2{\lm{z)\~^) . Therefore, 

aAr(z) = -TrRiv(z) + (66) 
where ei^Niz) is the holomorphic function on C\IR defined by 

e\,N[z) i CF ^ , „ cr,r-^ „ cr^ 



— Tr Ajv(z)Ra,(z) + -TrE [Qjv(z)] Rjv(z) — Tr AAr(z) 



|ei,jv('Z)l £ Pi(|z|)P2(|Im(z)r^). Finally, &n{z] canalso be writtenas 

a . eiiv(-z) 
aNiz) = -TrRiv(z) + (67) 
iV iv^ 

where ei^Niz) is equal to 

1 + CTajv(z) 

ei,Af(.z) = —e\,N{z), 

1 + (Ta;Af(z) 

and satisfies \ei^Niz)\ < Pi(|z|)P2(|Im(z)r^). 

In order to evaluate the behaviour of a^(z) - ■^TrR^(z), we differentiate (66) w.r.t. z and get the following 
result. 

Proposition 6. For z e CW, it holds that the derivatives e'^ (z) and e'-^ ^(z) ofe\^N{z) and ei,iv(z) w.r.t. z satisfy 
\e\ j^{z]\ < Pi(|z|)P2(|Im(z)ri) and\e[ j^(z)\ < Pi(|z|)P2(|Im(z)|-^). 

Proof: The proof uses Lemma [T4l and the observation that the spectral norms ||R^(z)|| and||R^(z)|| are bounded 
by Pi(|z|)P2(|Im(z)ri). The details are omitted. □ 
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In order to complete the proof of the lemma, we establish that 
Proposition 7. For z e CW, 



N2 ' 
N2 ' 



where \e2,Niz)\ and\£2,Niz)\ are both bounded by Pi {\z\]P2{\^^iz)\ 
Proof: We first observe that l|66) and (67] imply that 



'1,N 



iz) 



a'j^{z)-d'j^iz)^a-Ti(R'j^{z))-d'j^{z)+ „ , 
N N'^ 

1 - ^'l N^Z) 

a^(z) - 5^(z) = (7-Tr {B!f^{z)] - 8'f^{z) + 



We start with the classical identities 



Riv(z) -Tiv(z) = Rjv(z) (Tjv(z)-i -Rjv(z)-i)Tjv(z), 
Riv(z) - Tiv(z) = Riv(z) (f w(z)-i - RwCz)"!) j^^^^^ 



and get that 



a 



— Tr {Rjv(z) -TAf(z)) = (ajvCz) -6n{z)] zv^iz) + [aj^iz) -8n{z)) u^iz), 
■^Tr (Rjv(.z) - Tjv(z)) = (aAf(z) - 5iv(z)) Mjv(z) + [aN(z] -Sn{z))zvn(z], 



with 



and 



(68) 
(69) 



(70) 
(71) 



un{z) - — Tr 



Rw(z)BjvB;,Tjv(z) 
N " (1 + crajv(z))(l + o-<5iv(z))' 



Un(z] - — Tr 



RN{z)B*^BNtN{z) 



N (l + o-ajv(z))(l + o-5iv(z)) 



VNiz) = — TrRjv(z)TA,(z) 



5jv(z) = — TrRjv(z)Tjv(z). 



Note that it is easy to check that ujv(z) = un{z). We differentiate (70), (71)) w.r.t. z, we use (68) , (69) and Proposition 
m and recall that both |a]v(z) - 5]v(z)| and |ajv(z) - 5jv(z)| are bounded that i^Pi(|z|)P2 (|Im(z)ri) (see (23]). We 
check that u^iz), zv^iz), zv^iz) are their derivatives are bounded by Pi(|z|)P2 (|Im(z)r^), and obtain eventually 
that 



q:^(z)-5^(z) 
a^(z)-5^(z) 



Miv(z) ZL'jv(-Z) 

zi)iv(z) ujv(z) 



a^(z)-5^(z) 
a' (z)-5' (z) 



^2 



e3,iv(z) 
e3,iv(z) 



with |c3,jv(z)|, |e3,iv(z)| boundedby Pi(|z|)P2 (|Im(z)| We denote by A]v(z) the determinant of the above system, 
i.e. 



A]v(z) = (1- Uiv(z)) -zy]v(z)i)jv(z). 



(72) 



The determinant AAr(z) was studied in fT4l and in f23] where it was proved that |AAr(z) ^| < Pi(|z|)P2 (|Im(z)| ^) 
on a subset &n of C defined as 



@jv = |z e C - K, — Qi(|z|)Q2 (llm(z)r^) < 1 
where Qi and Q2 are 2 polynomials independent of N. Thus, we can invert the previous system on @]v to get 



a'j^iz) - d'j^iz) 


1 


l-Miv(z) ZVn{z) 


1 


es.Niz] 


a'^iz) - S'^lz) 


Ajv(z) 


ZVf^[{z] 1-M]v(z) 


N2 


£3,n{z] 
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This implies that |q;^(z) - 5'^(z)| is bounded by -p-PidzDPa (llmfz)!"!) on ®jv. If z £ C\{IRu®iv}, we use the trick 
in [12] . We remark that 

\a'^{z)-d'f^{z)\ < |a^(z)| + |5^(z)| < 

\lmz\ 

for each z, and that 1 < ^Qi (|z|)Q2 (|Im(z)r^) on C\ {Ru ®jv}. Therefore, 

\a'^{z)-5'^{z)\ < -^-Lqj(|z|)Q2 (llm(z)r^) < -^PidzDPa (|Im(z)ri) 
on C\ {R u ®Ar}. This in turn shows that (63) holds on C\R. □ 

5.2 Proof of lemmalS differentiability of ^Tr ^ (n n) 

We first need to establish the following useful Lemma. 

Lemma 16. Given an integer D > 0, let f be a continuous real function on U^. Let © he an open set ofR^ such 
that U^W has a zero Lebesgue measure. Assume that f is a^^ function on and that its gradient f on 6 can be 
continuously extended to K^. Then f is on the whole R^ with gradient f . 

Proof: We only need to prove that for any x £ R^ - ^ and any sequence x„ — ► x, 

fix,,)- fix) = {f'ix),Xn-x) + 0idn). 

where dn-\\Xn-x\\. Since / is uniformly continuous on any small neighborhood of x, there exists a sequence 5„ 
such that for every y and y' in this neighborhood for which ||y - y'|| < 6n, \ fiy) - fiy')\ ^ d^. Since -0 has a 
zero Lebesgue measure, there exists y„ and z„ in & such that 

lUn-ynll <mm{d n,dl) and ||x-z„|| <min(5„,d^). 

Therefore, it holds that max(|/U„)-/(y„)|,|/(z„)-/(x)|) < dl Writing /(xj-/(x) = /(x„)-/(y„)+/(y„)-/(z„) + 
fi^n) ~ fix), we obtain that - fix) - /(y„) - /(z„) + oid„). By differentiability of / on ^ and continuity of /' 
at X, 

fiyn)-fiZn) = {fiZn),yn- Zn) + oiWyn- ZnW) = </'(x), X„ - X> + 0((i„) 

which proves the lemma. □ 

We now complete the proof of the Lemma. We consider ^ £ Sf^(R,R), and establish that, considered as a 
function of the real and imaginary parts of Wjv, function j^Jr^' (CIm) is continuously differentiable on R^*'^^, i.e. 
that for each pair ii,j), the partial derivatives 

d (1 ~ , 

exist, and are continuous We denote by & the open subset of K^'^'^ for which the eigenvalues (A/ jv)z=i,...,m of 
ZjvZ^ have multiplicity 1. It is clear that M.'^'^'^W has a zero Lebesgue measure. On 0, it is standard that the 
eigenvalues (A;,jv)/=i,...,m are'i^^ functions and that 

au^-t^-"''-]/.'- '''' 

Using Lemma 4.6 in Haagerup-Thorbjornsen fT2l, we obtain 



dWi 



{JiVfiClN)} = Tr U'{n^)—^ — [CIm]] = 

i,;,JV I dWij^N I 



M 

z;vX[v>'("iv)]//n/,iv 

1=1 



(74) 



and get that ■pTr'I'(nAr) isa'^^ on 0. By Lemma [161 it remains to establish that the righthandside of (74) can be 
continuously extended to any point of R^*'^^\^. For this, we first prove the following useful result. 



^ is real valued, the partial derivatives w.r.t. Wi j pj thus coincide with the complex conjugate of the partial derivative w.r.t. Wj j pj. It is 
therefore sufficient to consider these derivatives. 
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Lemma 17. IfXk.N = M,N' then liff{ClN)]kk = [Vi^N)]ii- 

Proof: We start by observing that for any integers mi,m2,..., mt, matrix A = Aj^ 11 Ajy ■ ■ ■ 11 Ajy writes 



A = 





J mi 






J mi 


S mi 


"■M,N 





imt-i 



S mf-i 



1 m(-i 
^M,N 



:m, 



hence [A] fcjt = [A]// if A|t,jv = ^/,jv- The same can be said about ll^A and All^. Consequently, the result of the 
lemma is true when i/i" is a polynomial. Since any continuous function y) is the uniform limit of a sequence of 
polynomials on compact subsets of IR, the result is true for such -if/. □ 



We consider an element of R^Miv^ 



, and denote by mi, ... , mi, with M - Zflj mi, the respective multiplic- 
ities of the eigenvalues of where 2^ - Bjv + W^. We also denote by (11; 7v)/=i,...,l the orthogonal projection 
matrices over the corresponding eigenspaces. LemmafTTl implies that for each i - 1, . . . , L, 



mi+...+m, + m,-+i-l,mi+...+ mi+m,-+i-l 



J mi+...+ mi,mi+...+m,- 

Therefore, for any sequence (W^r ,i)„epy converging torward W?,, it holds that 



d 1 - , 
lim Tr>P njv] 



Wjv=Wjv,„ 



1=1 



J.I 



This completes the proof of Lemma[6l 

5.3 Proof of lemmaS uniform boundedness of E[||WAr||P] 

It is clear that it is sufficient to prove the boundedness of EiHWivP] if the entries of Wjv are real. We thus consider 
the case of real matrices and denote by the largest singular value of The following concentration result is 
well-known. 

Theorems (|7,Th. 11.13]). It holds that^[Xj^\ < 1 + and for all t> 0,P[Xn > l + ^/c^+ t] < exp[-Nf/2]. 
Using Theorem|6]and for p > 2 the inequality, 

n^oo n+oo 

E[X^]=j^ P>{XN>t)ptP-^dt<p{l + y/c^)'' + j^ P{XN>t+l + y/c^)plt+l + y/c^)P-'^dt, 
we easily obtain E[X^] < < oo, with K a constant independent of N, for all p e N. 

5.4 Proof of LemmaH} differentiability of the regularization factor 

We first establish that det 0(Zjv2;^) is a function, and that (50) holds. We use the same approach as in Haagerup 
& Thorbjornsen [12] Lem. 4.6] . We start begin by shov«ng that the differential of det cj){,X) is given by 

det 0(X)'.H = Tr (adj((/)(X))</)'(X)H) . (75) 

As det (X)'.H = Tr (adj(X)H) and (X")'.H = X'HX""i-' for any « e N, we have 

det (X")'.H = Tr (adj(X")(nX""i)H) 

since adj(X") and X commute. So (75) is true when is a polynomial. By choosing a sequence of polynomials P„ 
such that Pn ^ <p P'„ <p' uniformly on compact subsets of U, we generalize (75) to any cpe'^i. Now one can 
check that 



-,- e,e^-z.^, 



(76) 
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and it remains to apply the composition formula for differentials to obtain (50). 

We also remark that at a point Wjv for which there exists a A/ ]v t supp ((/)), we have 



^ i , ^ . 

adj(0(ZjvZjv*))(/)'(Zjv2:jv*) = ^ (]! ^Wfc.wjJ^'U/.jvju^u; = 



u 

:[ 

hence the derivative tSOl is zero on si^^^. 



It is easy to check that det (p{tljqtl^) is a 'i^i function on the open set @ of all matrices for which the 
eigenvalues of Zjv^^ are simple, and that (51) holds if Wjv e i.e. on a set of probability 1. In order to show 
that det (p[{lj\f^N) is a S^i function on we use again Lemma[T6] and verify that (51) can be continuously 

extended to tf'^'^W. For this, we claim that 

[adj((/)(niv))(/)'(njv)]j^j^= [adj((/)(niv))(/>'(«iv)],,; (77) 
if ^k,N = ^i,N- Indeed, given e > 0, let (j)e{x) - (^(x) + e. Since ^^(QAr) > 0, 

adj((^.^(njv))0U"iv) = det ((^^(njv))</);nnAf)(^';;(njv). 
Applying LemmaflTlto y/ - cj)'^^ ^ (t>'e' we obtain that 

[adj(0E(njv))(^^(nw)]fcfc= [adj((/),(nA,))(^;(njv)]„ if Afc 

and letting c ^ 0, we obtain the same result for adi{(p{ClN)](p' . Similarly to the proof of LemmalU this proves 
that (51) can be continuously extended to if-^^W. 

5.5 Proof of lemmadl} various estimates 

In this section, we denote by ar,iv(-z)>Sr,iv(z),Rr,jv('Z) and Rr,jv(-3) the regularized versions of the respective func- 
tions ajv(z),a]v(z),Riv(z) andRjv(z) defined in Section [5J] i.e. 

ttr.AfU) = crE|-^Tr(Q]v(z))X]vj and ar,iv(z) = (7E|-^Tr(Qjv(z))xjv 

and 

Rr,wU)= !L^-za + (Jar,Nizm ,Rr,N(z]^\- —-za + (Jar,Niz)]\ . 

\l + aar,Niz) I \l + aar,Niz) j 

It is clear that ar,N and fir.iv are the Stieltjes transforms of positive measures carried by C\supp((/)) and C* \supp((^) 
respectively and with mass gcnEIxn] and gEIxnI- This implies that the following uniform bounds hold: Let Jff 
and be compact subsets of C\supp((^) and C* \supp((^) respectively, then we have 

sup |Q;r,Jv(-z)l < C and sup |ar,Jv('Z)l < C (78) 

In order to establish Lemma[12l it is necessary to show that similar bounds hold for functions ^^^^ ^^^^ , ||R,-,]v(z) || 

and ||R,_jv(-z)ll. For this, we introduce function Wr,Niz) - z(l + (Tar,iv(-z))(l + crdtr^Niz)) and prove the following 
lemma 

Lemma 18. For any compact subset^ o/C\supp(0), it holds that 

sup |ar,iv(-z)-5jv(z)| '0, (79) 

zejr JV-oo 

inf min U^jv- "^r]v(-z)| > C> 0. (80) 

z^Xk=\ m' ' ' ' 
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Proof: Define )Cjv(z) ar,N(z) - Sn{z) where we recall that S^iz) - acMrriNiz) - j^Tr:{T]si{z)). Since S^iz) and 
ar,Ar(z) are Stieltjes transforms of positive measures carried by C\supp((/)), ?cjv is holomorphic on C\supp((/)) and 
satisfies 

C 

|7CJV(Z) 



d(z,supp((^)) 



This implies that the sequence {k^) is uniformly bounded on each compact subset of C\supp((^). By Montel's the- 
orem, (jcjv) is a normal family. Let (k:^(]v)) a subsequence of (?Cjv) which converges uniformly to k on each compact 
subset of C\supp ((/))) . Then k is holomorphic on C\supp((^). From |23 Prop.6], E [-^TrQ^vCz)] - -^TrTjvtz) — ► for 

z E CMR"*" and since xn ^jv 1 a.s., dominated convergence theorem implies 



7Cjv(z) = E [-^TrQjv(z)Xiv] - ^TrTiv(z) ^ 



-N '"J N 

for z e C\IR^. Thus, k(z) = for z e C\IR^, and by analytic continuation, jc(z) - for all z £ C\supp((^). Therefore, 
all converging subsequences extracted from the normal family [KNiz]] converge to uniformly on each compact 
subset of C\supp ((/)). Consequently, the whole sequence [kn) converges uniformly to on each compact subset of 
C\supp(0). This completes the proof of 179). We also notice that 

ail- Cm) a[l-C]v), , 

ar,N(z) = a,-,jv(z) - ^ — + — 1 - E(xjv) (81) 

z z 

and recall that 5]v(z) = SnIz] - ^lizM ^ _ ^^^^-^ ^ ^ g^^^-j^ integer p, (79) implies 

sup |z(ar,jv(z)-5jv(z))| ^0. 



Hence, it holds that 



sup I Wr,Niz] - Wjsiiz) \ 0. 



Thus, (80) follows immediately from (l3) . □ 

LemmafTsl immediately implies that the following uniform bounds hold. 

Lemmal9. LetJ^ andJ^ be compact subsets of C- supp{(p] andC* -supp{(p) respectively. PorN largeenough, we 
have 



sup 


1 


< c. 


(82) 


1 + aar,Niz) 




sup ||Rr,Jv(z)| 
Z£j*r 


<c. 


(83) 




sup ||Rr,Ar(z)| 

zejr 


<c. 


(84) 


sup||Rr,iv(z)-Tjv(z)| 


^0, 


(85) 


SUp||Rr,iv(z)-fjv(z)| 


^0. 


(86) 



Z£jff 

Proof: We first recall that inequality (10) holds. Therefore, the uniform convergence result (79) implies that 

inf \l + aarNiz) \ > - 
zejfr' ' ' 4 

for N large enough. This establishes (82) that holds for N large enough. In order to prove (83) , we express Rr ^v(z) 
as 

Rr,iv(z) = (l + crar,Jv(z))(BArB^- Wr,jv(z))" 
and use (78) and (80) . The proof of (84) is similar, and is based on the identity 

Rr,jv(z) = (l + (Tar,jv(z)) (B^Bjv - Wr,jv(z))~^ . 
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We remark that function &r,Niz) has a pole at z = 0. Hence, any compact J? over which ||Rr,jv(z) || is supposed to 
be uniformly bounded should not contain 0. The proof of (85) follows immediately from (79) and from (82) 
(84) . Finally, to establish (86) , we remark that 



Tr,jv(z) ■■ 



B^Rjv(z)B 



h 



Wr.Niz) l + (jar,Niz) 



l + cr5jv(z)' 



and that \ Wr,Niz)\ and Iwjvl'Z)! are uniformly bounded from below by (T3) and (80) (recall that is one of the 
eigenvalues of Ba/B^) . □ 

We now establish (56) and (57). In order to prove that ajvU) - (5jv(z) - &[^] on C\U^, |10| and |23| used 
the integration by parts formula (see e.g. [19]) and the Poincare inequality to show that the entries of E[Q]v(z)] 
are close from the entries of R]v(z) (see the fundamental equation (65)). Then, ajv(z) -5jv(z) was evaluated by 
solving a linear system whose determinant Ajv(z) given by (72) was shown to be bounded from below. Lemma [191 
allows to follow exactly the same approach to establish (56) and (57) . However, functions ajv, SAr,Rjv,R]v have to 
be replaced by their regularized versions. The following results show that the presence of the regularization term 
Xn does not modify essentially the calculations of |10| and |23|. We first indicate how the integration by parts 
formula is modified. Vec(.) denotes the column by column vectorization operator of a matrix. 

Lemma 20. Let (/jv) jv>i be a sequence of continuously differentiable functions defined on cM(M+iV) j^j-^/z polynomi- 
ally bounded partial derivatives satisfying the condition 

sup |/jv (Vec (Qjv(z)) , Vec(Ziv)) Xn\<C. 



Then, for all peN,we have 



E [/ (Vec (Qjv(z)) , Vec(i:iv)) liv] = E / (Vec (Qiv(z)) , Vec(Ziv)) x% + 



e \,N[z] 
NP 



for all keN* , and 



a2 






E 


N 










E 


N 





(5/(Vec(Qjv(z)),Vec(Ziv)) 



dW 



■Xn 



ij.N 



a/(Vec(QAf(z)),Vec( Xjv)) 



-Xn 



£2,Niz) 
NP ' 

£3,Niz) 

NP ' 



(87) 



(88) 



(89) 



E [ Wij,Nf (Vec (Qw(z)) , Vec(Ziv)) Xn\ = 

with sup^^g^^ |£-;,jv('Z)l £ C < oo. 

As for the use of the Poincare inequality, we have: 

Lemma 21. Let (Mjv(z)) a sequence of deterministic complex M x M matrix- valued functions defined on C\R such 
that 

sup ||MAr(z)|| < C. 



Then, 



sup Var 



TrQw(z)Mjv(z)ii, 



andforaN e C*^ such thatsup^ ||ajvll < cxd 



C 



C 



sup Var[a^Qjv(z)M]v(z)ajvXiv] < — . 



N 



Moreover, the same kind of uniform bounds still hold when Qjv(z) is replaced by Qjv(z)^. 
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The proofs of these results are based on elementary arguments, and are thus omitted. Following the calcula- 
tions of [10] and |23], we obtain that 



E [QNiz)XN] = Rr,jv(z) + A,,jv(z)Rr,jv(z) + E [Qjv(z)ljv] R,,iv(z)— Tr Ar,jv(z) + 0A,(z)Rr,iv(z) (90) 

for each z e C\supp((^) where 0Af(z) is a matrix whose elements are uniformly bounded on dSSy by for each p, 
and where Ar,w(z) is the regularized version of matrix Ajv(z) introduced in lemma[T5]defined by 



Ar,jv(-Z) 



1 



rE[Qjv(z)iJv]E 



— TrQjv(z)xjv-E 



(J 

— TrQjv(z)Xiv 



— TrZ^Q]v(z)BivXiv 



1 



1 + crar,jv(z) 
1 



(J 



— TrZ^Qiv(z)BivX]v-E 



— TrZ^Qiv(z)BjvXjv 



N 



^TrQjv(z)X]v-E 



(T 

— TrQjv(z)Xiv 



l + aar,N{z) 

After some calculations using Lemmas [19] |20l|2T] we eventually obtain that 

C 



Qjv(z)Z]v2:^Xiv 



(91) 



sup |a^(E[Q]v(z)xjv]-Rr,iv(2))aiv| 



^3/2' 



I o" I C 

sup |ar,]v(.z)- — Tr(Rr,iv(2:))| < 

I a . \ C 

sup |a;r,w(z)-— Tr(Rr,iv(z))| < 



ZEl3S?„ 



(92) 
(93) 



for all large N. In order to prove (56) and (57), it remains to handle the terms involving the difference Rr,iv(z) 
T]v(z). We show in the following that 



C 



sup \a*^ (Rr,]v(z) - T]v(z)) ajv| < 



(94) 



for all large N. We start as usual with the identity Rr,Ar(z) - Tjv(z) = Rr,N(z) (T]v(z)"i - Rr,Ar(z)"i) x^(z), to get 

a,-jv(.z)-5jv(z) 



(Rr,Af (z) - Tjv(z)) aAf = cj 



at,Rr,w(2:)BjvBt,TAf(z)aiv 



(1 + aar,N{z)) (1 + cr(5iv(z)) "~ 
+ ZCT (ar,Ar(z) - <5]v(z)) a^R^.Af (z)T]v(z)ajv. 

The expression (M) of ar,N implies that z(q:,-,jv(z) - Sm{z)) - ziar,Niz) - 5iv(z)) + {-^) for each integer p. Thus, 
to prove (56) and (57) , it is sufficient to check that 



sup |ar,]v(z)-<5jv(z)| 



c_ 

N2- 



We will use the same ideas as in Section [Ol and remark that (a,-,jv(-z) - Sn{z}, &r,Nlz) - 5jv(z)) can be interpreted 
as the solution of a 2 x 2 linear system whose determinant is a regularized version of (72) , and appears uniformly 
bounded away from zero on dS^y. 

Using again the previous expression of Rr,]v(z) - Tjv(z) together with (92) . (93) and repeating the procedure for 
Rr.Niz) -TAr(z), we obtain 



ar,N(z)-SNiz) 
ar,Niz)-§Niz) 



Ur,Niz) ZVr,Niz) 
ZVrNiz) UrNiz) 



ar,N{z)-8N{z) 
ar,Niz)-SN{z) 



1 



cjv(z) 
£nIz) 



(95) 



withM,,jv(z) = ^Tr^j^ 



Rr,jv(z)BLBjvTjv(z) 



,-]v(z))(i+g5jv(z)) ' ^r.Niz) = ^TrRr,iv(z)TAr(z) and £'r,jv(z) = ^TrRr,jv(z)TAf(z). The quanti- 
ties ejv(z), ejv(z) satisfy sup^^g^ kivt'Z)! < C,sup^^g^ |cjv(z)| < C. The determinant of the system is given by 



Ar.Af (Z) = (l - Ur,N {z)f - Vr.Af (z) Vr,Niz) . 
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Lemma [T9l implies that for all large N, Ur,Nlz), Vr,Nlz) and Dr,Niz) are uniformly bounded on dS^y. Therefore, to 
conclude the proof of (94), it remains to check that for aU large N, 

inf |Ar,iv(z)| > C>0. 

Consider the function Ajv(z) where we have replaced the matrix Rr.jvtz) and R,,jv(-z) byTjv(z) and Tjv(-z),i.e 

Ajv(z) = (1 - UNiz)f - vm{z) DNiz), 



with iijv(z) = ^Tr 



jv [1 Jffi^'°H(i+a5w(z)) ' '^nIz) = ^TrTiv(z)2, and vn(z) = ^TiTNizf- Denote by hN(z) = An(z) - 
An(z). Lemmas [Tsl [191 imply that |ujv(z)- ujv(z)|, \vn{z]- vn{z]\ and |i)]v(z)- 5jv(z)| converge to uniformly on 
9SSy which of course implies 



sup |/ijv(z)| 



N- 



0. 



(96) 



Using Cauchy-Schwarz inequality, we get 

|Ajv(z)| > A^(z) (1 - u^(z))^ -\z\^Vj^[z)v^[z), 

with M^(z) = ^Tr^^^^^^^^gj^, y^(z) = ^TrTiv(z)Tiv(z)* and £^(z) = ^Trf iv(z)Tiv(z)*. Now, we use the 



HN^^' - N 

following lemma. 



Lemma 22. There exists a constant C > independent ofN such that 

inf A^(z)>C. 



Proof: It is shown in 1 23 1 and 1 14 1 that A^ (z) is the determinant of the following 2x2 linear system 



Im(5jv(z)) 
Im (z5jv(z)) 



Uj^iz) Vj^iz) 

\z\^D^[z) u^{z] 



Im (5jv(z)) 
Im (z5jv(z)) 



Im(z) 



Rn'^z) 
u^[z) 



(97) 



and that for z £ 
easily get that 



, Ajy(z) > 0. Solving the system, and looking at the corresponding expression of Im(5jv(z)), we 



Im(z) a 
^iv(^) = T , ,, T:TrTAf(z)Tiv(z)*. 
Im(5jv(z))Af 



for z e C\IR. Expressing Tjv(z) - Tiv(z)* as TAr(z) (Tjv(z)"* -Tjv(z)"i)Tjv(z)*, and using the equation 5jv(z) = 



§Tr(T]v(z)), we obtain that 



and deduce the useful formula 



Im(5jv(z)) = Im(Wiv(z)) -TrTjv(z)Tjv(z)' 
N 



4iv(^) ■ 



Im(z) 



Im(WAr(z)) 



(98) 



Using the integral representation (5iv(z) - J^^ — 
pression i^'w^jz)) extends to C\5^Ar, and therefore to C\supp((^), and satisfies 



we obtain after straightforward computations that the ex- 



|Im(wjv(z))| 

sup < C 

z€d<^, |Im(z)| 



Eq. (98) thus implies that 



sup |A^(z)| <C. 



which concludes the proof. 



□ 
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We deduce from this that inf^^^^^ | Ajv(z)| > C > for all large N. Therefore, we can invert the system (95) and 
obtain 

C 

sup \ar,Niz)-SM{z)\< 

Z£d3ty 

for all large N. This establishes (57) and completes the proof of (56). 

The proof of (58) is similar to the proof of LemmaH) but as above, aNiz),aN{z),RNiz) and Riv(z) have to be 
replaced by their regularized versions ar,Niz), ar,Niz),Rr,Niz) and Rr,NN{z). The reader can check that the prop- 
erties of these regularized functions aUow to follow the various steps of the proof of LemmaH) 
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